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ABSTRACT 

Faint submillimeter sources detected with the Submillimeter Common-User Bolometer Array on the 
James Clerk Maxwell Telescope have faced an identification problem due to the telescope's broad beam 
profile. Here we propose a new method to identify such submillimeter sources with a mid-infrared image 
having a finer point spread function. The Infrared Space Observatory has provided a very deep 6.7 //m 
image of the Hawaii Deep Field SSA13. All three faint 850 /zm sources in this field have their 6.7 /xm 
counterparts. They have been identified with interacting galaxy pairs in optical images. These pairs 
are also detected in the radio. Two of them are optically faint and very red (I > 24, I — K > 4), 
one of which has a hard X-ray detection with the Chandra satellite. As these observing properties are 
similar to those of local ultraluminous infrared galaxies, their photometric redshifts are derived based 
on submillimeter to mid-infrared flux ratios assuming a spectral energy distribution (SED) of Arp220. 
Other photometric redshifts are obtained via \ 2 minimization between the available photometry data and 
template SEDs. Both estimates are in the range z = 1-2, in good agreement with a spectroscopic redshift 
and a millimetric one. The reconstructed Arp220 SEDs with these redshift estimates are consistent with 
all the photometry data except Chandra's hard X-ray detection. The sources would be a few times 
more luminous than Arp220. With an assumption that AGN contributions are negligible, it appears 
that extremely high star formation rates are occurring in galaxies at high redshifts with massive stellar 
contents already in place. 
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1. INTRODUCTION 



The installation of the Submillimeter Common-User 
Bolometer Array (SCUBA) on the 15 m James Clerk 
Maxwell Telescope has brought us the discovery of a new 
population of faint submillimeter sources (Smail, Ivison, & 
Blain 1997; Barger et al. 1998; Hughes et al. 1998; Eales et 
al. 1999). Meanwhile, the far-infrared background (FIRB) 
was detected with the Far Infrared Absolute Spectrom- 
eter (FIRAS) and the Diffuse Infrared Background Ex- 
periment (DIRBE) on the COBE satellite (Puget et al. 
1996; Schlegel, Finkbeiner, & Davis 1998; Hauser et al. 
1998; Fixsen et al. 1998). The FIRB is comparable to or 
larger than that in the optical wavelengths, and much of it 
could be explained by integrating individual submillime- 
ter sources down to the faintest limit (Blain et al. 1999). 
Thus, it is indispensable to investigate the faint submil- 
limeter population in order to understand overall history 
of energy production in the Universe. 

Although SCUBA has unveiled the faint submillimeter 
population, its beam is too broad (15" FWHM at 850 /xm) 
to pinpoint their optical counterparts (e.g. Barger, Cowie, 
& Richards 2000). The high density of the submillimeter 
sources also raises a problem of source confusion (Eales et 
al. 2000; Hogg 2001). Up until now, methods to bypass the 



broad submillimeter beam have included centimetric ra- 
dio interferometry with VLA (Ivison et al. 1998; Richards 
1999; Smail et al. 1999; Barger et al. 2000) and millimetric 
interferometry with IRAM or OVRO (Downes et al. 1999; 
Gear et al. 2000; Bertoldi et al. 2000; Frayer et al. 2000). 
Most recently, an attempt to observe at a longer wave- 
length has started at the IRAM 30 m telescope (Bertoldi 
et al. 2000). The larger dish provides a slightly smaller 
beam (10" at 1.3mm). 

Here we investigate the identification of SCUBA sources 
with a mid-infrared deep image having a relatively small 
beam (7'.'2 at 6.7 /im). Faint submillimeter sources, pre- 
sumably dusty systems at high redshifts, could be bright in 
the mid-infrared regardless of their energy sources. If the 
submillimeter emission was produced by dust heated by 
star forming activity, there should be stellar systems fol- 
lowing the stellar initial mass function. The stellar systems 
could be luminous at the rest-frame near- infrared (> 1 fim) 
even in the presence of dust. In the case that the source 
was powered by an active galactic nucleus (AGN), very 
hot dust surrounding the AGN could emit its reprocessed 
energy at > 2 /zm. For both cases, the emitted light could 
be received at the mid-infrared, even though the sources 
were at high redshifts. 



1 Based on observations with ISO, an ESA project with instruments funded by ESA Member States (especially the PI countries: France, 
Germany, the Netherlands and the United Kingdom) and with the participation of ISAS and NASA. 

2 Institute of Astronomy, University of Tokyo, 2—21-1 Osawa, Mitaka, Tokyo, 181-0015 Japan; ysatoSioa. s .u-tokyo . ac . jp 

3 Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI 96822 

4 Astronomical Institute, Graduate School of Science, Tohoku University, Aramaki, Aoba, Sendai, 980—8578 Japan 

5 Institute of Space and Astronautical Science (ISAS), 3-1-1 Yoshinodai, Sagamihara, Kanagawa, 229-8510 Japan 

6 Gunma Astronomical Observatory, 6860-86 Nakayama, Takayama, Agatsuma, Gunma, 377-0702 Japan 

1 



2 



Y. Sato et al. 



2. SOURCE IDENTIFICATION 

A sensitive mid-infrared camera ISOCAM (Cesarsky et 
al. 1996) on board the Infrared Space Observatory (ISO, 
Kessler et al. 1996) provided us with the first opportu- 
nity to identify faint submillimeter sources at the mid- 
infrared. In the Hawaii Deep Field SSA13, a very deep 
mid-infrared map was obtained with the LW2 (6.7 /im) 
filter (Sato et al. 2001). This field was also surveyed ex- 
tensively with SCUBA (Barger et al. 1998; Barger, Cowie, 
& Sanders 1999a). In the upper left panel of Fig. 1, we 
plot both 6.7/xm and 850 /xm sources detected within the 
SCUBA field-of-view. All three 850 /im sources; A, B, 
and C (these identifiers are used throughout this paper; 
Table 1) are found to be located very close to certain 
6.7 /xm sources. Monte Carlo simulations of source de- 
tection suggest that displacement of the detected position 
as large as the beam size could happen in low signal-to- 
noise ratio (S/N) ranges (Eales et al. 2000; Hogg 2001). 
We estimate the probability of random alignments be- 
tween the 6.7 /im and 850 /xm sources. As the surface den- 
sity of 6.7 /xm sources at this depth is larger than that 
of 850 /im sources, we assume a Poisson distribution for 
6.7/xm sources. Then, the chance alignment probabil- 
ity can be defined as P = cxp(— 7rr 2 7V(6.7/tm)), where 
r is the distance between 6.7/xm and 850 /tm sources and 
iV(6.7/xm) is a 6.7/xm integral number count at the flux 
of the corresponding 6.7 /im source. The chance probabil- 
ity became only P — 0.9, 16, and 4.1 % for sources A, B, 
and C, respectively. Then, we associate the closest 6.7 /im 
sources with the corresponding submillimeter sources. 

The optical counterparts of the sources were identified 
in deep HST images taken with the WFPC2 F814W filter 
(Cowie, Hu, & Songaila 1995). Interacting galaxy pairs 
were found at the peaks of the 6.7/xm contours for all the 
850 /tm sources (Fig. 1). The occurrence of an irregular 
morphology is consistent with the results of ultraluminous 
infrared galaxies (ULIGs) discovered with IRAS (Sanders 
& Mirabel 1996). Within 1" from all the optical counter- 
parts, faint 1.4 GHz sources are detected (E. A. Richards 
2001, private communication). Source B also has a hard 
X-ray (2-10 keV) detection with Chandra (Mushotzky et 
al. 2000). The small rate of AGNs in submillimeter sources 
(1/3) is similar to the previous results (Almaini, Lawrence, 
& Boyle 1999; Fabian et al. 2000; Hornschemeier et al. 
2000; Severgnini et al. 2000; Barger et al. 2001). Magni- 
tudes at the K, I, and B bands were measured with ad- 
ditional data taken from the ground (Cowie et al. 1996). 
Two of the three sources are optically faint and very red 
(Table 1). Here are notes for the individual sources. A) 
This is the most significant SCUBA detection in this sam- 
ple (4.7 a). The ISOCAM contour reduces the uncertainty 
of the SCUBA position. In the middle of the 6.7/xm con- 
tour, there is a faint (/ ~ 25) interacting galaxy pair hav- 
ing a very red color of / — K > 4. B) One of the 6.7 /im 
peaks matches a colliding galaxy pair with / ~ 24 and 
I — K > A. The ISOCAM source appears somewhat ex- 
tended probably due to confusion of faint sources. C) A 
bright merger (/ ~ 22) with tidal tails is located at the 
edge of the beam centered at the nominal SCUBA posi- 
tion. The original SCUBA image shows a sign of elonga- 
tion, suggesting a lower position accuracy (Fig. 1 in Barger 
et al. 1998). 



3. REDSHIFT ESTIMATES 

The three SCUBA sources detected at 6.7/xm share the 
same observational properties as those of local ULIGs; i.e. 
an irregular morphology, a red color, and optical and X- 
ray faintness. Therefore it seems possible to estimate their 
rcdshifts with an assumption that their spectral energy 
distributions (SEDs) are similar to that of Arp220, an 
archetypal ULIG in the local universe. We adopt a UV 
to submillimeter SED in the GRASIL library to represent 
Arp220 (Silva et al. 1998). Applying a convolution with 
two bandpass filters at 850 /im and 6.7 /xm, we derived flux 
ratios /„(850 /tm) / /„(6.7 /xm) as a function of redshift (up- 
per left panel of Fig. 2). Because of a monotonic increase 
with redshift, this relation can be used to estimate red- 
shifts, analogous to the redshift estimator based on a sub- 
millimeter to radio flux ratio (Carilli & Yun 1999, 2000; 
Barger et al. 2000). The derived estimates were in the 
range z p h(6.7/xm) = 1-2 (Table 2). The reconstructed 
SEDs are consistent with other photometry data, suggest- 
ing the usefulness of this redshift estimator (Fig. 2). 

With the identification results above, we could omit the 
assumption of an Arp220 SED to obtain photometric red- 
shifts. Other estimates z p h(SED) were derived via \ 2 min- 
imization between all the photometry data and SED tem- 
plates. We here adopt the GRASIL SED library (Silva et 
al. 1998); UV to submillimeter SEDs of local templates 
(starbursts: M82, NGC6090, and Arp220, normal galax- 
ies: M51, M100, and NGC6946, and a giant elliptical) and 
UV to radio SEDs of evolving E, Sa, Sb, and Sc galaxies. 
We add an X-ray to radio SED of a heavily obscured AGN 
NGC6240 compiled by Hasinger (2000). The best scaling 
factor for a certain SED at each redshift was determined 
by minimizing x 2 . The \ 2 value was summed up at the 
wavelengths where the model predictions could be evalu- 
ated. Upper limits were used as restrictions. Thus, the 
number of degrees of freedom is the number of detections 
(within the wavelength range of the SED) minus one (for a 
scaling factor). The smallest reduced \ 2 was achieved with 
the Arp220 SED for sources A and B, and with the 13 Gyr 
Sc SED for source C. The resulting z p h(SED)-values are 
similar to z p h (6.7/xm). A spectroscopic redshift of source 
C was measured as z sp = 1.038 (Cowie et al. 1996). A 
millimetric redshift of z p h(L4 GHz) = 1.8 was derived for 
source B (Barger et al. 2001). All of these are consistent in 
the range z — 1-2 (Table 2). Source C gives a minimum \ 2 
at z ~ 1 for many SEDs including the Arp220 SED. The 
hard X-ray flux for the Chandra source B was not used in 
the x 2 calculations except in the case with the NGC6240 
SED. An overall fit to the hard X-ray to radio data turned 
out to not be so good; the NGC6240 SED suggests z = 1.9 
but with a reduced % 2 — 5.9. (cf. Fig. 11 in Barger et al. 
2001). If we admit such a level of reduced x 2 5 there could 
be some room for a higher redshift for the optically faint 
sources A and B. Source A could be at z = 5.3 with the 
0.8 Gyr E SED (a reduced X 2 = 1-9). For source B, the 
0.2 Gyr E SED gives z = 3.1 with a reduced X 2 = 3.2 (for 
the UV to radio data only). Both E SEDs represent a 
dusty elliptical before the galactic wind blows. 

4. DISCUSSION 

We identified source C with a galaxy offsetting 6" 4 
from the nominal SCUBA position. Such sizes of offsets 
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were observed in simulations of the source detection in a 
crowded field (Eales et al. 2000; Hogg 2001). SCUBA ob- 
servations of radio selected sources actually found such 
offsets at the submillimctcr flux level of 2-4 mJy with 
S/N = 3-5 (Barger et al. 2000). The expected large offsets 
and optical faintness suggest that the mid-infrared identifi- 
cation of SCUBA sources could be an effective way to find 
their optical counterparts. For the submillimeter sample 
presented here, their counterparts are as faint as 14^iJy at 
6.7/xm, K = 20, or I = 25. Integral galaxy counts at the 
corresponding depths are 1 x 10 4 deg -2 , 3 x 10 4 deg -2 , 
and 8 x 10 4 deg -2 , respectively (Sato et al. 2001; Maihara 
et al. 2001; Metcalfe et al. 2001). Thus, the probability of 
finding spurious associations by chance would be smaller 
by a factor of 3-8 at 6.7/zm. 

Our ISOCAM map taken with a pixel field of view 
(PFOV) of 6" still has a low resolution problem. The 
6.7 /zm contour for source C extends to the east to in- 
clude an edge-on spiral. The extended appearance is also 
seen for source B. Although their fluxes were corrected 
to be one point source statistically (Sato et al. 2001), ef- 
fects of the extension to the corrected flux for a particular 
source remain uncertain. A mid-infrared imaging with a 
finer PFOV (l"-2") with SIRTF/IRAC or ASTRO-F/IRC 
would clarify such ambiguities. The high sensitivity of 
these detectors will also strengthen redshift estimates for 
the milli-Jansky 850 /im sources. They can be detected 
even at z = 5, and their redshifts could be determined, 
provided that the Arp220 SED is a representative for them. 
For cases of other SEDs, multi-band mid-infrared fluxes 
should be required to estimate photometric redshifts be- 
cause such high redshift sources are likely to be faint in 
the optical, or even in the near-infrared. 

The estimated range of redshift for the three sources 
in the SSA13 field is z = 1-2. Given a cosmology with 
Hq = 50kms -1 Mpc -1 and Qq = 0.1, the sources would 
be 2-4 times more luminous than Arp220. These prop- 
erties are similar to other submillimeter sources (Lilly et 
al. 1999; Barger et al. 1999b, 2000). A simple scaling of 
the model parameters for the GRASIL Arp220 SED gives 
star formation rates in the range (1-2) x 10 3 M Q yr -1 and 
stellar masses in the range (3-7) x 10 11 M Q . While the 



peak of the dust emission was not observed directly, the 
6.7 /im flux originated from the rest-frame near-infrared. 
Because rest-frame near-infrared emission has a good cor- 
relation with stellar mass, it would be accurate to assume 
that these submillimeter sources have massive stellar pop- 
ulations in them. With signs of a major merger in all the 
optical counterparts, similar to the irregular morphologies 
already reported (Smail et al. 1998; Lilly et al. 1999), it 
is suggested that less massive galaxies identified at z ~ 3 
(Pcttini et al. 1998) would evolve into local massive galax- 
ies via major merging with the submillimeter bright phase. 

5. CONCLUSIONS 

All three of the faint 850 /im sources (2-4 mJy) in 
the SCUBA SSA13 deep field have been found to have 
ISOCAM 6.7 ^m counterparts (10-30 /zJy). Utilizing the 
smaller beam size at 6.7 ^m, we find that all the three 
sources are coincident with interacting galaxy pairs in the 
HST I band images. They all have VLA 1.4 GHz coun- 
terparts, and only one of them is detected with Chandra 
at the 2-10 keV band. Based on the properties similar to 
local ULIGs, we used submillimeter to mid-infrared flux 
ratios to estimate photometric redshifts. We derived other 
photometric redshifts utilizing the X-ray to radio photom- 
etry data. Both estimates are consistent with a spectro- 
scopic redshift and a millimetric one. An Arp220 SED at 
z = 1-2 provides a good fit to the available data. With an 
assumption that AGN contributions are negligible, the im- 
plied SFRs arc in the range (12) x 10 3 M yr -1 and the 
mid-infrared emission requires stellar masses in the range 
(3-7) x 10 n M o . 
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Table 1 



Identified submillimeter sources 



ID 


Name 1 


RA 


Dec 


2-10 keV 2 


B 




K 


6.7 /im 


850 ^tm 1 


1.4 GHz 4 




[BCS99] 


(J2000) 


(J2000) 




(mag) 


(mag) 


(mag) 


(MJy) 


(mjy) 




A 


J131232. 1+424430 


13 12 31.94 


+42 44 29.7 


<3.2 


>25.5 


>24.7 


19.7 


18 


3.8 




B 


J131228.0+424458 


13 12 28.29 


+42 44 54.6 


5.67 


>25.5 


24.0 


19.9 


14 


2.3 


39 


C 


J131225. 7+424350 


13 12 25.18 


+42 43 44.9 


<3.2 


24.3 


21.6 


18.5 


27 


2.4 





References. — (1) Barger et al. 1999a; (2) Mushotzky et al. 2000; (3) Sato et al. 2001; (4) Barger et al. 2001 



Note. — Right Ascension in hour, minute, second and declination in degree, arcmin, arcsec are measured at the K band. 
Hard X-ray fluxes are in units of 10 ergs cm . For the magnitudes, 3" diameter aperture measurements are corrected 
with average offsets to approximate total magnitudes (Cowie et al. 1996). Upper and lower limits are 3<r. 



MIR identification of faint SMM sources 



5 




J AJ VJ ■ 7 T. J k t , "Tat, ■ 



Fig. 1. — Identification of submillimeter sources A, B, and C (IDs in Table 1). In the upper left panel, the spatial distribution of the 
submillimeter and mid-infrared sources is shown in the SCUBA SSA13 deep field (2.'7 in diameter; dotted circle). North is up and east to 
the left. SCUBA 850 (im and ISOCAM 6.7 (im sources (Barger et al. 1999a; Sato ct al. 2001) arc marked with large and small circles with 
their beam sizes (15" and 7"2 FWHM, respectively). Magnified views of the submillimeter positions are displayed on the iTST/WFPC2 / 
band images (20" width in gray scale). The 6.7 um circles are substituted with contours to make full use of the smaller beam size and higher 
significance of the 6.7 fim sources than those of the 850 ^tm detections. Contour levels are (2.5, 3.0, and 3.5) a, while (2.0, 2.5, and 3.0) a for 
source B. Three VLA 1.4 GHz sources and one Chandra 2— lOkeV source are shown with crosses and a diamond, respectively (E. A. Richards 
2001, private communication; Mushotzky et al. 2000). 
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Table 2 

Redshifts of the submillimeter sources 



ID V z p h(SED) z ph (6.7^m) z ph (1.4GHz) 2 



A 


2 +0 - 4 


2 +0 - 6 

z - u -0.6 




B 


i-e±8:i 


i -°-0.6 


1.8 


C 1.038 


o.9±g:i 


i o +0 - 3 

i - u -0.3 





References. — (1) Cowie et al. 1996; (2) Barger et al. 2001 

Note. ■ Photometric redshifts z p h(SED) are derived via 
X 2 minimization between the available photometry data and 
template SEDs. The uncertainty ranges are 1 a confidential 
limits with the best fit SED. Other estimates z p h(6.7/zm) are 
obtained with the Arp220 SED (Fig. 2). The la errors origi- 
nate from 1 a photometry errors at 6.7 and 850 fim. 




-4 -2 2 4 6 -4 -2 2 4 6 

Log (Wavelength [|im]) Log (Wavelength [|im]) 



Fig. 2. — Submillimeter to mid-infrared flux ratios /„(850 £tm)//„ (6.7 fim) for an Arp220 SED, showing a monotonic increase with redshift 
(upper left panel). Based on this relation, photometric redshifts 2 p h(6.7 (im) are derived for the three submillimeter sources A, B, and C. 
With these redshift estimates, Arp220 SEDs are reconstructed to fit the 850 fim and 6.7 fim measurements. We adopt the GRASIL SED (Silva 
et al. 1998), complemented with radio and X-ray data (Carilli & Yun 2000; Iwasawa 1999). Reconstructed SEDs arc consistent with other 
photometry data. The hard X-ray detection of source B is beyond the wavelength span of the SED. The fluxes and upper limits (3 a) are 
measurements at 2-10 kcV, 0.5-2 keV, B, I, K, 6.7//m, 450 fim, 850 /urn, and 1.4 GHz (Table 1; Mushotzky et al. 2000; Barger et al. 1998). 



